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Abstract 

A periodically forced oscillator in a model for seasonality shows stability pockets and chains thereof 
in the parameter plane. The frequency of the oscillator and the season indicated by a value between 
zero and one are the two parameters. The present study is intended as a theoretical complement to 
the numerical study of Schmal et al. [5] of stability pockets or Arnol’d onions in their terminology. 

We construct the Poincare map of the forced oscillator and show that the Arnol’d tongues are taken 
into stability pockets by a map with a number of folds. Stability pockets are already observed in an 
article by van der Pol & Strutt in 1928, see and later explained by Broer & Levi in 1995, see [3]. 

Keywords: forced oscillator, resonance, synchronization, Poincare map, stability pocket, circadian clock, 
Zeitgeber, entrainment 


1 Introduction 

The numerical study by Schmal et al. of a model for seasonal effects on the circadian clock shows 
stability pockets and chains thereof in the parameter plane. This study aims to complement their results 
by giving a theoretical background for the observed phenomena. We also indicate what is to be expected 
when their system is perturbed. The setting of the problem is bifurcations of parameter dependent 
dynamical systems, in particular periodically forced oscillators. 

A periodically forced oscillator shows periodic dynamics or synchronization if the frequency of the forcing 
is close enough to a rational multiple of the frequency of the oscillator in absence of forcing. On the 
other hand, if this frequency ratio is not close enough to a rational number for many frquency ratio’s 
the dynamics is quasi periodic. Now we consider the periodically forced oscillator as a system of two 
asymmetrically coupled oscillators, where the first is forced by the second. The latter has fixed dynamics. 
In particular its frequency is fixed and we set it equal to one by a scaling of time. Thus the frequency 
of the first oscillator becomes an essential parameter of the system. Another important parameter is the 
coupling strength. In the parameter plane of coupling strength versus frequency there are regions, called 
Arnol’d tongues (see m) or resonance regions, where the first oscillator synchronizes with the second, 
that is there are one or more stable (and unstable) periodic solutions. These tongues are wedge shaped 
and their vertices lie at rational points on the frequency axis. So for a fixed value of the coupling strength 
there is a frequency interval where synchronization occurs. On the boundary of this interval the periodic 
solutions disappear in saddle-node bifurcations. If the coupling strength goes to zero the interval shrinks 
to a (rational) point. If the coupling strength increases there may be many other bifurcations, for example 
period doubling bifurcations, see |4j. 

In certain examples of periodically forced oscillators the tongues close again in a second vertex, forming 
a so called called stability pocket (called Arnol’d onion by Schmal et al. [5]), see left part of figure]^ 
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This phenomenon occurs for example in Hill’s equation and it has already been observed by van der Pol 
& Strutt [7], although they do not mention it explicitly as such. Much later Hill’s equation has been 
reanalyzed by Broer & Levi [3] using methods not available to van der Pol & Strutt. Then the term 
instability pocket was introduced. However, it depends on the point of view whether a pocket is called a 
stability or an instability pocket, see remark [3731 Although the context of the Hill equation differs from 
ours, the mechanism by which the pockets are formed is the same. In both systems a general underlying 
system has a wedge shaped resonance region. The map that takes the parameter space of the specific 
example to the parameter space of the general system has one or more folds, thus generating pockets. 

The model in Schmal et al. [5] consists of a differential equation for a specific two dimensional oscillator, 
in fact the normal form of the Hopf bifurcation, with periodic forcing. Here we take a slightly different 
approach. Our oscillator will be a general phase oscillator asymmetrically coupled to a second phase os¬ 
cillator with constant frequency, the forcing oscillator or Zeitgeber. First we study a general coupling and 
approximate the Poincare map. Then we use these results to study seasonality like in by considering 
a periodic block function as forcing, where the length of the block (modeling the length of the daylight 
interval and thus season) is controlled by an additional parameter. 

We could have started with a general oscillator in two dimensions with an external periodic forcing. But 
then the technicalities of bringing the system in manageable form and keeping it so when the external 
forcing is applied would obscure the phenomenon we want to study. Our primary goal is the way the 
dynamics in the phase direction of the oscillator changes when the external force is applied: we are less 
interested in the change of the shape and position of the periodic orbit. Therefore we start with a phase 
oscillator from the very beginning. 


2 Two asymmetrically coupled phase oscillators 


Let 4’ Etnd ip be phase angles so '0; ‘z? G <5'^ where ip is interpreted as the phase angle of the oscillator 
and ip as the phase angle of the external forcing. Furthermore, let be a smooth function on T^. The 
following asymmetrically coupled system describes a periodically forced phase oscillator. 

{4 = uj + pf^[4,ip) 

I 0 = 1 


The system depends on small parameters pL,5, - ■ ■ G M™. The frequency of the oscillator is uj ior p — 0. In 
that case the dynamics of the system is quasi-periodic if u! is irrational, but if u is rational, the dynamics 
is periodic. Thus for /i = 0 periodic dynamics is exceptional and quasi periodic dynamics is typical. This 
changes dramatically if /r 7 ^ 0, then periodic dynamics becomes typical. Thus it seems natural to take a 
closer look near rational values of w, therefore we set ui = | with p and q relative prime integers and 
(5 is a small parameter. It is most efficient to study this system via a Poincare map. A good candidate is 
the map that scores tp at consecutive crossings of (/? = 0, sometimes called the stroboscopic map. Since 
all Poincare maps are equivalent, results will not depend on this choice. For convenience we switch to 
the lift of this system to then the differential equation becomes 

(x=^+6 + fif^ix,y) 

1 y = 1 

Now we define the (lift of the) Poincare map as follows. 

Definition 2.1. Let <I>t be the Bow of equation then the Poincare map F is implicitly defined by 
{F{x),l) = ^4x,0). 


The map F in this definition is the lift of a circle map / (not to be confused with /^), the Poincare map 
on the circle. In actual computations it is easier to use the lift F. Now since F is the lift of a circle map 
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it must be of the form F{x) = x + uj + php^{x) where /i is a 1-periodic function. Unfortunately there is 
no easy connection between h and / in equation ([^. But by a so called normal form transformation we 
obtain a vectorfield approximation of the Poincare map F, see section]^ The approximating vectorfield 
is in the righthand side of the following equation 

« = + (3) 

where fk,i,m depends in a complicated way on the Taylor-Fourier coefficients of and the sum runs over 
all k and I with pk + ql = 0 and m S {1, • • • , n}. To be more precise, the time-1 flow of this vector held 
approximates the Poincare map. This is formulated more precisely in section [4. 1| 

The main property of the Poincare map / we use is that a stationary or a periodic point of / corresponds to 
a periodic solution of equation Q and vice verse. Thus the existence of periodic solutions of equation 0 
can be read off from the existence of stationary or periodic points of /. The vectorheld approximation does 
not cover the full dynamics of the original system. But hyperbolic stationary points of (|^ correspond to 
stationary points of the Poincare map F. Moreover by persistence of saddle-node bifurcations we recover 
the well-known picture of resonance tongues for F in the (w,^)-plane from equation ([^. For each pair 
{p,q) a tongue emanates from the w-axis at w = U The tongue at {p,q) = (1,1) is called the main 
tongue, the others are just labeled by the pair {p, q). Indeed we may rewrite the differential equation as 
follows 

u = 5 + pff,{u) (4) 

where is a 1-periodic function. Then at least near (0,0) stationary points exist in a region in the 
(i 5 , /x)-plane bounded by curves of saddle-node bifurcations of the form p = 71 (5 and p = 72(5, for some 
constants 71 and 72. In the simplest case where f^{u) = sin(27ru) we have p = ±( 5 . 

In the next section we take a specific form for equation ([^ containing parameters a and A and we study 
the map (ct. A) 1 —> {5, p). The inverse image of this map takes the tongues of the general equation for two 
asymmetrically coupled oscillators to the resonance regions of model with season dependent forcing. 


3 Main results 

3.1 Stability pockets in the seasonal oscillator 

In a simple model of season dependent forcing we take a phase oscillator with a particular 1-periodic 
forcing. We will call it the seasonal oscillator for short. The general form will be the following. 

ii: = uj + r]fi'ilj)+egx{ip) 

\(p = l 

Here / and gx are C°° functions on so that the corresponding vector field is also . The function / 
describes the non-linearity of the oscillator and gx determines the external periodic forcing. As before we 
take (jj = ^ + a, with positive integers p and q. The parameters a, 77 and e are small, but not necessarily 
of the same order. 

To model seasonality we use a function as in namely one that depends on a non-small parameter A, 
which determines the fraction of the period the forcing is ’on’. Thus we interpret A = 0 as ’winter’ and 
A = 1 as ’summer’. A simple example being a block function with a block of length A S [0,1] 




1, 0 < t < A 
0, A < t < 1 


( 6 ) 
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with g\{t + 1) = g\{t) for all t. This function is not C°°. Therefore we repl ace g\ by the convolution 
(j> * g\ where (j) is for example the function (p^x) = y^exp(—In section 4.3 we will show that the 
result does not depend on the choice of (p. For simplicity we assume from now on that g\ is C°°. 

We use the methods of the previous section to find a vectorfield approximation of the Poincare map for 
the seasonal oscillator. The following theorem gives a more precise statement. For a proof see section [4^ 
where also the meaning of degree will be clarified. 

Theorem 3.1. The vectorfield approximation up to degree two of the Poincare map of equation ([^ is 


u = a — 




(7) 


where u = x — The Fourier coefhcients of / and g\ are fk and gkW respectively, where 

i 

5 fc(A) = exp(-—) —(exp(-27rifcA) - l). 


We wish to find a region in the parameter space of the seasonal oscillator, that is in the (a, A)-plane, 
where stationary points of equation Q exist. Since this equation is a special form of equation ([^ we 
consider the map (cr, A) i-P- (S,fj,). The result is formulated in the following theorem. 

Theorem 3.2. Let p and q be relative prime as mentioned before. Then the map (cr, A) i—?> {S,p) is 
implicitly given by: 

{5, fj.) = {a- g^ci - er]Xc 2 , er]\h{X)\). 

where ci and C 2 are constants and h is a ^-periodic function with zero average. This implies that the 
(p, q)-tongue has p stability pockets. 


The map implicitly defined in the theorem takes the main tongue of a general periodically forced oscillator 
into a so called stability pocket in the (a, A)-plane of the seasonal oscillator. Other tongues, depending 
on p, are mapped to a chain of stability pockets. In figure [T] a graphical representation of this result is 
shown. 








Figure 1: Stability pockets and folds of the map (cr. A) i—>■ {S,fj,). Left: the pocket for the {l,l)-tongue; 
right: the chain of pockets for the {2,l)-tongue. 
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Remark 3 . 3 . 


i) Theorems 3.1 and 3.2 together support the numerical results in [S], there is an excellent qualitative 
agreement with their findings and the propositions. There is no reason to doubt that with some 
more effort this agreement can also be made quantitative. 


ii) In [5] the stability pocket is called “Arnol’d onion”. Here the term stability pocket is used to connect 
with the existing literature on Hill’s equation, see [3]. In Hill’s equation the zero solution becomes 
unstable in the resonance tongues. In case such a tongue closes, a “pocket” is formed which is then 
called an instability pocket. In our case we are interested in synchronization, that is in the existence 
of stable periodic solutions. Since these exist in resonance tongues, a tongue closing and forming a 
“pocket” is now called a stability pocket. 

iii) In the special case = sin(27r'0) the map becomes simpler, namely {S,p) = (cr — rj^ci — 

sr]Xc2,er|^/2\ sin(p7rA)|). In fact the map (d, p) = {a — X,\ sin(p7rA)|) has essentially the same proper¬ 
ties. This map is used to draw the graphs in figure 

iv) We have to find the map (cr. A) i—>■ {S, p) for each tongue. The reason is that no single parameterization 
of all normal forms exists: for each resonance p : q we have to compute a new normal form. 


v) In a more accurate model of seasonality, the parameter A will vary in time. Since this variation 
is slow with respect to the frequency of the forcing we may consider a model based on a slow-fast 
system. From this perspective we are studying the fast limit of such a system. 


In a biological interpretation of these results, the image of the main tongue is perhaps the most relevant, 
see figure For a fixed value of A at a certain season, the range of entrainment is indicated by a 
horizontal line intersecting the stability pocket. Only in this interval of uj values, the seasonal oscillator 
entraines to the Zeitgeber. That is for these w values, the oscillator has periodic solutions with the same 
period (or frequency) as the Zeitgeber. From the perspective of the oscillator the value of oj is fixed. Then 
the range of entrainment is in the seasonal direction. As can be seen in the figure, the oscillator cannot 
entrain to the Zeitgeber for all seasons. In case of the main tongue, not in the middle of summer and 
not in the middle of winter. When we look at the (2,1)-tongue, the range of entrainment in the seasonal 
direction may consist of two intervals, see figure For more biological interpretation and examples refer 
to [5]. 


1 


A 

0 




2 u) 


Figure 2: Stability pockets and ranges of entrainment in both frequency and seasonal directions. The 
image of the main tongue is shown in the left figure. The right figure shows the image of the (2, l)-tongue 
for a small perturbation of g\. 


Proof of theorem \3.S\ Compare equation @ to equation 0 to determine the map {5, p) (cr, A). We 
immediately see that d = a — rfci — eriXc 2 with ci = p fkf-k and C 2 = |/o. Recall that go = X. 
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The second component is determined by /mqS-mp exp(27rim(jfu). The A dependent factor of the 
modulus of the Fourier coefficients is | sin(TOp7rA)|. The latter is a periodic function of A with p + 1 
zeroes for all m. Or, put differently, with at least p folds for all m. Thus the number of stability 
pockets is determined by the first and largest Fourier coefficients, /i and /_i, and thus by the map 
{S, p) = (cr — — e? 7 Ac 2 , ejycal sin(p7rA)|). Therefore the (p, ( 7 )-tongues have a chain of p stability 

pockets. To describe their shape in detail one would need the remaining Fourier coefficients. □ 

3.2 Generalizations 

Let us now generalize the forcing gx to a function h\ such that still hg = 0 and hi = 1, so that A = 0 
corresponds to winter and A = 1 to summer. But for A € (0,1), hx is not necessarily a block function. 
For simplicity we assume that hx is C°°. Just like the Fourier coefficients of gx, those of hx can be 
regarded as parameterizations of a closed curve in C with endpoints at zero, the parameter being A. Now 
for fc 7 ^ 0 the A dependent factor of Pfc(A) is (exp(—27riA:A) — l), which passes through zero for A = ^ 
for n G {0,1,..., fc}. When we perturb gx to hx subject to the condition that hg = 0 and hi = 1, 
the Fourier coefficients again parameterize closed curves in C with endpoints at zero, but in general not 
passing through zero for A S (0,1). Thus we will still have stability pockets for the more general system, 
but the chains of stability pockets will open up to a single pocket, see the image of the (2, l)-tongue in 
figure 


4 Proofs 

This section contains the proofs of the theorems in the previous sections. We first give a short review of 
approximating a Poincare map and then apply this to our coupled oscillators. In a last section we discuss 
how to deal with a nonsmooth forcing. 

4.1 A vectorfield approximation of the Poincare map 

In this section we consider vectorfields on the phase space as lifts of vectorfields on the two torus. 
Our aim is to construct a vectorfield approximation of the Poincare map. Where the latter is the lift of 
the Poincare map on the two torus. We will work in the C°° context unless explicitly stated otherwise. 

First we define a parameter dependent differential equation on Let be a function —)■ M, 

1-periodic in both arguments and assume / has a Taylor-Fourier expansion 

Ux,y)= Y, 

k,l,m 

with coefficients Now consider the differential equation 

(i=^+S + pf^{x,y) 

\ y = 1 

depending on parameters S and fi. The first will be interpreted as a detuning and the second as the 
strength of the nonlinearity. The Poincare map F of this system is defined as {F{x), I) = $1(3::, 0 ), where 
<I>t is the flow over time t of equation Q. 

Computing the vectorfield approximation. By a sequence of near identity transformations we 
put system Q in a form such that the y dependence in the first component becomes trivial, at least 
up to a certain order. Here we use perturbation theory so we have to be precise about the size of 
various components. The grading we use will come from the parameters only. The reason is that the 
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vectorfield is periodic in x and y therefore we set degree(a:) = degree(?/) = 0. Furthermore we set 
degree(/i) = degree(i5) = 1. With these definitions the function / has a formal Fourier-Taylor expansion 
in X, y and p 

Mx,y)= 

k,l,m 

with Fourier-Taylor coefficients fk,i,m- Finally we define the vectorfield 


^0 = ^^ + ^ 


X = Xo + Xi + - 

^ _9 

^ dx~^ dy 
Xj = fk,l,r. 


+ Xn + ■ 


^m^2-!Ti(kx+ly) 


dx 


Now we apply a standard normalization procedure to obtain a normal form of the vectorfield X. The 
result is formulated in the next proposition. 

Proposition 4.1. By a sequence of n near-identity transformations the vectorfield X is transformed into 
X with 


X = Xo + Xi + --- + Xn + 0{\pr+^) 

\m\—j,pk-\-ql—0 

The new coefficients fk,i,m depend on the original coefficients in a complicated way. One exception being 
Xi for which fk.i^m = fk,i,m- The transformed vectorfield up to order n only contains resonant terms, 
that is in general fk,i,m ^ 0 only if the index (k, I, m) satisfies the resonance condition pk + ql = 0 . 

Proof. This will only be a sketch of the proof, for more details see [2]- The key idea of the proof is that 
every (7°° near identity coordinate transformation can be approximated as closely as needed by the flow 
of a C°° vectorfield [S]. Since we work in the context of C°° vectorfields on the two torus (or lifts thereof) 
the transformation has to respect this property. The flow of a C°° vectorfield on the two torus is such 
a transformation. Since these vectorfields form a Lie algebra we can use rather standard normal form 
theory. Our aim is to get rid of the time (y) dependence without transforming time, therefore we apply 
an asymmetric transformation generated by the vectorfield Y = g{x, y)-§^ so that the new x depends on y 
but not the other way around. The procedure is inductive by degree. Let X be the set of C°° vectorfields 
on the two torus depending on one or more small parameters p. We define a grading by the degree of p: 
Xk <Z X is the set of vectorfields of degree k or higher in p. In the normal form procedure we frequently 
use the commutator [•,•] of vectorhelds, then for X^ € Xm we have [Xm,X„] G Xm+n- Furthermore if 
Ym G Xm is a fixed vectorheld and adl)n : X [Ym,X] then {adYm)^{Xn) G Xkms-n- Because of these 
relations we may normalize for increasing degree. 

Now suppose we have normalized up to degree n — 1, then take Y of degree n. The transformation acts 
on the vectorfield X as expadT, namely 

expadT(W) = = Xq Xi -I- • • • -I- X^—i + X^ + adTT h.o.t. 

where adF(a;) = \Y, X]. As usual in normal form theory we try to solve Xn-\-adY (Aq) = 0 for g. Let gk,i be 
the Fourier coefficients of g then on the level of coefficients the equation becomes qfk,i—2'Ki{jpk-\-ql)gk,i = 0. 
Thus we set gk,i = 2 TTilpk'+qi) Provided that pk ql 7 ^ 0. This implies that in the normal form the so 
called resonant terms with coefficient fk,i satisfying the resonance condition pk-\-ql = 0 are retained. □ 
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Assuming we have normalized the vectorfield up to sufficiently high order, we only keep terms up to order 
n by truncation. As a last step we again change coordinates: u = x — and v = y. Then the new 
vectorfield Z becomes 

j=l \m\=i,pk+ql=0 

Indeed, \ipk+ql = 0 then changing to coordinates u and v we get kx+ly = k(u+^y)+ly = = 

ku. 

Interpretation of vectorfield Z. Now let X be the truncated normal form of X and Z be the 
vectorfield defined above. The corresponding flows are denoted by respectively $, $ and 4'. Let F 
be the Poincare map of $ defined as (F(x),l) = <i>i(x,0). In a similar way we define F, then this F 
is an approximation up to order n of F. Furthermore let ip be the flow of the first component of Z. 
Then we have 4's(u,0) = {ips{u),s) but also Ts(u,0) = (na;$s(u,0) — s^,s), where Ha, is the projection 

Tlx{x, y) = X. This means that F^(u) = n 2 ;$g(u, 0) = ipq[u) +p. Thus F‘i — p is equal to the time q flow 
of the first component of Z. Therefore we call Z a vectorfield approximation of F. 

Using the vectorfield approximation. The vectorfield Z approximates the Poincare map of the 
vectorfield A in Q but there is no equivalence. Therefore we have to be careful when drawing conclusions 
about the dynamics of (|^ from analysis of the vectorfield Z. Here we will be mainly interested in 
stationary points of Z. A first observation is that hyperbolic stationary points of Z correspond to 
relative equilibria of X and thus to periodic orbits of X (provided that the difference between X and X 
is small enough). Stationary points of Z satisfy equation 

71 

^ = ^ + E (9) 

i=l \m\=j,pk+ql=0 


Given p and q, the right hand side is a ^ periodic function. Therefore solutions, if they exist, come in q 
pairs. Now suppose solutions exist, then upon varying parameters p they may disappear in tangencies. 
Assuming for simplicity a single parameter p, we find a wedge shaped region in parameter space dehned 
by 7 i/i < (5 < 72 /i where solutions exist. This inequality only holds near {S,p) = (0,0) and the constants 
7i and 72 depend on the Taylor-Fourier coefficients fk,i,m- Dynamically speaking the boundaries of 
the wedge are curves of saddle-node bifurcations. These form the familiar stability tongues or Arnol’d 
tongues. Since saddle-node bifurcations in one parameter families persist under small perturbations, the 
stability tongues of vectorfield Z approximate those of vectorfield X. 

For each combination of p and q (positive and relative prime) we find a tongue. The one for {p,q) = (1,1) 
is called the main tongue. If we assume for simplicity that the Taylor-Fourier coefficients rapidly decrease 
for increasing k and I then the leading terms in equation 0 are 

0 = S + yp sin(27ru -|- y) (10) 


where 7 and y are determined by and /_iyy. This little calculation at least shows that the tip 

of the tongue is a straight cone. 


Similarly the tongues for {p,q) = {l,q) are determined by an equation whose leading terms are as in 
equation (10), but now 7 and y are determined by fq,-i,i and /_qy^i. Again this shows that the tip of 
the tongue is a straight cone, but the angle of the cone decreases with increasing q. 


Remark 4.2. We collect some remarks about the vectorfield approximation. 


i) The transformation u = x — ^y and v = y is not just useful but originates from the following. By 

construction the truncated vectorfield X = Xq + Xi -|- • • • -|- A„ commutes with Xq implying that 
their flows $ and also commute. In other words the flow of Xq generates a symmetry group 
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of X. By switching to ’co-moving’ coordinates, thus by the transformation the vectorfield X 
transforms to X — Xq. Therefore stationary points of X — Xq correspond to periodic orbits of X. 
Because of this property such stationary points are called relative equilibria. 

ii) In the normal form of the seasonal oscillator we will only need terms up to degree two. Then we 
have to compute a few higher order (h.o.t.) terms. Let X and Y be as in the proof, but now Y is of 
degree 1 then the terms we need are 

expady(X) = Xq + X^ + ady(Xo) + ady(Xi) + ^{adY)^{Xo). 


4.2 Approximation of Poincare map of the seasonal oscillator 


We apply proposition |4.1| to the differential equation of the seasonal oscillator. The starting point is the 
general form of the lift of equation (HI 

j X = u + r]f{x) + egx{y) 

\ y = 1 

We normalize this vectorfield to obtain an approximate Poincare map. We assume that / and g both 
have a Fourier series with coefficients fk and gk. Then the vectorfield is 


X = Xo- 
X - P ^ - 

q ox 


X1+X2 
dy 


Xi = {gf^x)+eg{y))^ = ^^ 


dx 


X2 = 0 — 
ax 

where fk,o = rifk and foj = egi. For both fc 7 ^ 0 and Z 7 ^ 0 we set fk.i = 0. In view of the resonance 


condition pk - 


-ql = 0 for a term exp(27rz(A:a; + Zj/))see proposition 


4.1 


of degree 1 , therefore we set degree(5) = 2 and degree(£:) = degree( 77 ) = 


there will be no resonant terms 
.. Note that cv = ^ + 6. 

The proof of theorem 3.1 consists of determining the vectorfield of degree 1 and computing terms up to 
degree 2 . 


Proof. Following the proof of proposition 4.1 we take a vectorfield Y of degree 1 and we try to solve 
Xi +ady(Xo) = 0 for the Fourier coefficients of Y. The degree is determined by the parameters <5, e and 
g. Let Y = a{x,y)-^ and a has Fourier coefficients ak,i. Then on the level of Fourier coefficients we get 
o-kfl = Vpy^fk if fc 7 ^ 0 and Uq,; = if ^ 0- Since we have solved equation Xi + adF(X q) = 0 we 

have ady(Xi) + f(ady)^(Xo) = —^ady(Xi). Thus we obtain up to degree two 

X = expady(X) = Xq + Xi + adX (Xq) + ady(Xi) + i(ady)^(Xo) + h.o.t. 


= Xq + 2 ^^■(" h-O-f■ 


As a shorthand write Xi = f-^,Y = ag| and ady(Xi) = h-^ then h = —/|y. From this expression 

we select the resonant terms, that is the hk^i satisfying pk + ql = 0 , then we are left with 


P P 

k pk-\-ql—0 

^ k m 
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Setting u = X — and v = y the vectorfield approximation of the Poincare map becomes 

< ^ k m 

1^ n = 1 

□ 


With this result we find the stability tongues of the vectorfield approximation. If {p, q) = (1,1) then the 
tongue boundaries follow from solving equation 


0 = (5 - ^ fkf-k f-<!k9ke 


27riqku 


for u, see section 4.1 


4.3 The Fourier coefficients of the forcing 

The forcing g\ of the seasonal oscillator, see equation is a piecewise constant function and therefore 
not C°° as required in proposition |4.1| on the approximation of the Poincare map. In order to get a 
smooth approximation of gx we use convolution with a so called Schwartz function. For this smooth 
approximation we find the Fourier coefficients. We are in particular interested in the dependence on the 
parameter A. The proofs of the following proposition and lemmas are found by straightforward arguments 
and computations. 

Proposition 4.3. Let (j> be a normalized Schwartz function and let (p be its Fourier transform. Then 
the Fourier coefficients of the convolution <f> * gx are p(k) ■ gk where go = 0 a,nd gk = ^ ^g-27rifeA _ ^j jf 
k ^ 0. In particular we have \gk \ = sin(7rA:A)|. 

Let ^ be a normalized Schwartz function, that is 

i) for all positive integers m and n we have sup^ < oo, 

ii) lR(l^ix)dx = 1. 

The familiar Gauss function (j){x) = ^e““^ is an example of such a function. For all a > 0 the function 
4>a{x) = ap^ax) is again a normalized Schwartz function and the larger a, the closer (pa * gx is to gx. 

Lemma 4.4. Let g be a possibly non-smooth 1-periodic function with a Fourier series such that g{x) = 
gk expipinkx) and let p be a normalized Schwartz function with Fourier transform p. Then p * g is a 
smooth 1-periodic function which has a Fourier series with coefficients p{k) ■ gk. 

This lemma shows that the results in proposition |3.2| do not depend on the choice of the Schwartz function 
p since the zeroes of (p{k) ■ 5 fc(A), as a function of A, are exactly those of 5 fc(A). 

By an elementary calculation we immediately hnd the Fourier coefficients of gx and a smooth approxi¬ 
mation of gx- 

Lemma 4.5. The Fourier coefficients gk of the function gx as defined in equation ([^, are go = X and 
~ k ^ 0. The Fourier transform of pa is pa{y) = exp(— 

Our main interest is in the A dependence of the Fourier coefhcients. So far we found p{k) ■ gk = 
exp(—^^^^)^(exp(—27rz/cA) — l) if A: ^ 0, from which we infer that 

\4’{k) ■ gk\ = e“'^^;^A/2(I - cos(27rfcA)) = e“^ sin(7rfcA)|. 
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